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Abstract 

Characterizing and quantifying quantum correlations in states of many-particle systems is at the core 
of a full understanding of phase transitions in matter. In this work, we continue our investigation of the 
notion of generalized entanglement [Barnum et al, Phys. Rev. A 68, 032308 (2003)] by focusing on a 
simple Lie-algebraic measure of purity of a quantum state relative to an observable set. For the algebra 
of local observables on multi-qubit systems, the resulting local purity measure is equivalent to a recently 
introduced global entanglement measure [Meyer and Wallach, J. Math. Phys. 43, 4273 (2002)]. In the 
condensed-matter setting, the notion of Lie-algebraic purity is exploited to identify and characterize the 
quantum phase transitions present in two exactly solvable models: the Lipkin-Meshkov-Glick model, and 
the spin- 1/2 anisotropic XY model in a transverse magnetic field. For the latter, we argue that a natural 
fermionic observable-set arising after the Jordan-Wigner transformation, better characterizes the transition 
than alternative measures based on qubits. This illustrates the usefulness of going beyond the standard 
subsystem-based framework while providing a global disorder parameter for this model. Our results show 
how generalized entanglement leads to useful tools for distinguishing between the ordered and disordered 
phases in the case of broken symmetry quantum phase transitions. Additional implications and possible 
extensions of concepts to other systems of interest in condensed matter physics are also discussed. 
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I. INTRODUCTION 



Quantum Phase Transitions (QPTs) are qualitative changes occuring in the properties of the 
ground state of a many-body system due to modifications either in the interactions among its 
constituents or in their interactions with an external probe |Q], while the system remains at zero 
temperature. Typically, such changes are induced as a parameter g in the system Hamiltonian H(g) 
is varied across a point at which the transition is made from one quantum phase to a different 
one. Often some correlation length diverges at this point, in which case the latter is called a 
quantum critical point. Because thermal fluctuations are inhibited, QPTs are purely driven by 
quantum fluctuations: fluctuations or correlations in the value of some observable or observables 
that occur in a pure state. Thus, these are purely quantum phenomena: a classical system in a pure 
state cannot exhibit correlations. Prominent examples of QPTs are the quantum paramagnet to 
ferromagnet transition occurring in Ising spin systems under an external transverse magnetic field 
J2H2H4L the superconductor to insulator transition in high-T c superconducting systems, and the 
superfluid to Mott insulator transition originally predicted for liquid helium and recently observed 
in ultracold atomic gases 

Since entanglement is a property inherent to quantum states and intimately related to quantum 
correlations |g], one would expect that, in some appropriately defined sense, the entanglement 
present in the ground state undergoes a substantial change across a point where a QPT occurs. Re- 
cently, several authors attempted to better understand QPTs by studying the behavior of different 
measures of entanglement in the ground state of exactly solvable models (see 
for representative contributions). Such investigations primarily focused on characteri zing entan- 
glement using information-theoretic concepts, such as the entropy of entanglement II 1 311 or the 
concurrence [14], developed for bipartite systems. In particular, a detailed analysis of the two- 
spin concurrence has been carried out for the XY model in a transverse magnetic field |Q, Q], 
whereas the entanglement between a block of nearby spins and the rest of the chain has been con- 
sidered in iflfll . While a variety of suggestive results emerge from such studies, in general a full 
characterization of the quantum correlations near and at a quantum critical point will not be pos- 
sible solely in terms of bipartite entanglement. Identifying the entanglement measure or measures 
that best capture the relevant properties close to criticality, including the critical exponents and 
universality class of the transition, remains open problems in quantum information and condensed 
matter theory. 
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InRefs. OQ], we introduced Generalized Entanglement (GE) as a notion extending the 
essential properties of entanglement beyond the conventional subsystem-based framework. This 
notion is general in the sense that it is definable relative to any distinguished subset ofobservables, 
without explicit reference to subsystems, which makes it directly applicable to any algebraic lan- 
guage used to describe the system (fermions, bosons, spins, etc.) Founding the 
notion on a distinguished set of observables makes it especially well suited to studying QPTs, as 
our definition makes the existence of GE equivalent to the existence of nonzero correlations or 
fluctuations in those observables. The basic idea is that any quantum state gives rise to a reduced 
state on the distinguished subset of observables [19]. These reduced states form a convex set; as 



with standard quantum states, there are pure (extremal) and mixed (non-extremal) ones [20]. We 
define a generalized entangled pure state, relative to a subspace of observables, to be one whose 
reduced state on that subspace is mixed. Although we will have little occasion in the present con- 
text to apply it to states that are mixed relative to the full set of observables, we extend this notion 
to include mixed states by defining a generalized entangled mixed state to be one that cannot be 
written as a convex combination of generalized-unentangled pure states. 

The special case in which the observable set is a Lie algebra is often important in Physics. In 
a broad class of such algebras described below, the algebra is not only a subspace of operators, 
but is such that we can define a natural Hermitian projection onto that subspace. Then a simple 
(global) measure of GE for quantum states is provided by what we call the purity relative to the 
algebra. This is defined as the squared length of the projection of the Hermitian operator (density 
matrix) representing the state onto the algebra. As argued in [1 1], if the correct algebra is chosen, 
the purity contains information about the relevant quantum correlations that uniquely identify and 
characterize QPTs of the system. 

In this paper, we deepen and expand the analysis initiated in Jill , by focusing on the detection 
of QPTs due to a broken symmetry as revealed by the behavior of an appropriate relative purity of 
the ground state. In Section|IIl we recall the relevant mathematical setting and the definition of the 
relative purity as a function of the expectation values of the distinguished observables. In Section 
Inil we discuss several examples where the relative purity is seen to provide a natural measure of 
entanglement. In Section |lVl we illustrate some physical criteria that are relevant in choosing the 
appropriate observable subalgebra and using GE as an indicator of QPTs. In Sections IVl and fVTl 
we explicitly characterize the QPTs present in the so-called Lipkin-Meshkov-Glick (LMG) model 
I21H22I and in the one-dimensional spin- 1/2 anisotropic XY model in a transverse magnetic field, 
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respectively. This is done by studying the properties of the purity relative to different algebras 
of observables in the ground state of both models. We find the relevant critical exponents for 
these models, and in the case of the anisotropic XY model in a transverse magnetic field, obtain 
a new "global" disorder parameter, the variance of the number of spinless fermionic excitations 
in a Jordan- Wigner-transformed representation of the system. Finally, we provide in separate 
Appendices the details underlying various statements made in the main body of the paper. These 
include the relationship between standard separability and GE (Appendix A), the GE properties 
of two special classes of spin states, the cluster and valence bond solid states (Appendix B), the 
proof of the relationship between the local purity and the Meyer- Wallach entanglement measure 
(Appendix C), and the semiclassical properties of the LMG model in the thermodynamic limit 
(Appendix D). 



II. GENERALIZED ENTANGLEMENT AND RELATIVE PURITY 



In the GE approach, entanglement is considered as an observer-dependent property of a quan- 
tum state, which is determined by the physically relevant point of view through the expectation 
values of a distinguished subset of observables. Whenever a preferred decomposition into subsys- 



tems is specified in terms of an appropriate (physical or encoded ||23 



24, 



25 



2611 ) tensor prod- 



uct structure, GE becomes identical to standard entanglement provided that distinguished observ- 
ables corresponding to all local actions on the individual subsystems are chosen: in particular, for 
H = ®{Hi with dimCHi) = di, standard entanglement of states in H is recovered as GE relative to 
tyioc — ®iSu(di) JTli ll5ri (see also Appendix lA"b. In fact, the subsystems relative to which standard 
entanglement is defined (whether directly identifiable with physical degrees of freedom or related 
to "encoded" or "virtual" ones) are always understandable in terms o f approp riate (associative) 



algebras of local observables. This has been observed before, e.g. in ||25 



2611 (see also 12711 for 



a recent analysis). However, it is important to realize that the GE notion genuinely extends the 
standard entanglement definition, and does not coincide with or reduce to it in general. On one 
hand, this may be appreciated by noticing that even for situations where a subsystem partition is 
naturally present, states which are manifestly separable relative to such a partition may possess GE 
relative to an algebra different from h; oc (see the two spin-1 example of Section lITlh. On the other 
hand, as also emphasized in Jlli . GE is operationally meaningful in situations where no phys- 
ically accessible decomposition into subsystems exists, thus making conventional entanglement 
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not directly definable. 



A. Relative purity for faithfully represented Lie algebras 

As mentioned in the Introduction, we will focus on the case where the distinguished observables 
form a Lie algebra f) of linear operators, acting on a finite-dimensional state space H for the 
system of interest, S. (Note that we will not usually distinguish between the abstract Lie algebra 
isomorphic to f), and the concrete Lie algebra f} of operators that faithfully represents it on H.) We 
will assume h to be a real Lie algebra consisting of Hermitian operators, with the bracket of two 
linear operators X and Y being given by 

[X,Y]=i{XY-YX). (1) 

In this way, operators in () can be directly associated with physical observables. For the same 
reason, we will also use a slightly nonstandard (but familiar to physicists) notion of the Lie group 
generated by h, involving the map X e lX instead of the mathematicians' X i— > e x , for X e 
h. No assumption that the Lie algebra acts irreducibly on H (i.e., that it admits no nontrivial 
invariant subspaces) will be made, but important consequences of making such an assumption will 
be discussed. We will also assume the Lie algebra to be closed under Hermitian conjugation. This 
implies that it is a reductive algebra (not to be confused with reducibility of the representation). 
In our context, a reductive Lie algebra is best thought of as the product (direct sum, as a vector 
space) of a finite number of simple Lie algebras, and a finite number of copies of a one-dimensional 
Abelian Lie algebra. A simple Lie algebra is a non-Abelian one possessing no nontrivial ideals, 
where an ideal is a subalgebra invariant under commutation with anything in the algebra; the 
relevant property of ideals here is that they can be quotiented out of the algebra, allowing it to be 
written as a nontrivial product of ideal and quotient; thus simple Lie algebras are non-Abelian ones 
that cannot be decomposed into factors, so the factorization used in defining reductive Lie algebras 
above is maximal. The product (direct vector-space sum) of a finite number of simple Lie algebras 
is called semisimple, and thus a reductive algebra is the product of a semisimple and an Abelian 
part. The reader is referred to lEsl 2^1 I30I I31I for relevant background on Lie algebras and their 
representation theory. As this subsection unfolds, we will summarize much of this representation 
theory in a way suited to our needs, and the reader should concentrate on understanding the content 
of the statements, and not vex him or herself unduly about understanding why they are true. 
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We will consider pure quantum states of S, \ip) E H, as well as mixed quantum states of 
S, described by density matrices p acting on H. Since f) was assumed closed under Hermitian 
conjugation, the projection of a quantum state p onto f) with respect to the trace inner product 
is uniquely defined. Let denote the projection map, p i— > V^(p). If p is a pure state, p = 
\ip) {ip | , the purity of\ifj) relative to f) (or t) -purity) is defined as the squared length of the projection 



according to the trace inner product norm 111 51 : that is 



W» = Tr[(^>^|)) a ]. (2) 

The f)-purity may be explicitly evaluated upon selecting an operator basis B = {A±, . . . , A L } for 
f). By assuming the A a to be Hermitian, 



and orthogonal, 



A a = Al, (3) 



Jr(A a Ap) = 5 a ,p , (4) 



Eq. © may be rewritten as 

L L 

= Tr [ E MA aP )MAp P )A a Ap] = J2( A *) 2 > (5) 

a,P=l a=l 

where (A a ) denotes the expectation value of the observable A a in the pure state j-0). 

An important property following is that the f)-purity is invariant under group transformations: 

L 

if a new basis for f) is introduced by letting A a = D^A a D, with D = exp(i tpAp), D^D = 1, 

(3=1 

and tp real numbers, then one finds 

L L 

hw) = E^) 2 = E^-) 2 = • (6) 

a=l a=l 

Sometimes it is useful to introduce a common normalization factor K in order to set the maximum 
value of the purity to 1, in which case Eq. © becomes 

L 

JW» = »<5>«> a . (7) 

a=l 

As mentioned earlier, a pure quantum state [ip] is defined to be generalized entangled (general- 
ized unentangled) relative to (} if it induces a mixed (pure) state on that set of observables. When 
P) is a complex semisimple Lie algebra acting irreducibly on H, it was shown in I15I (Theorem 



14, part (4)) that \ip) is generalized unentangled with respect to t) if and only if it has maximum 
h-purity, and generalized entangled otherwise. Under the same assumptions, the abovementioned 
Theorem (part (3)) also leads to the identification of the generalized unentangled pure states as the 
generalized coherent states (GCSs) associated with h |32, In other words, all generalized 

unentangled states are in the (unique) orbit of a minimum weight state of f) (taken as a reference 
state) under the action of the Lie group. Remarkably, GCSs are known to possess minimum in- 



[35 



36D, so that, similar to the familiar 



variant uncertainty, (AF) 2 (\ip)) = J2 a i^a) ~ (A a ) 2 
harmonic-oscillator ones, they may be regarded in some sense as closest to "classical" states. 

Our characterization theorem for generalized unentangled states on irreducible representations 
used some standard facts from the theory of semisimple Lie algebras and their representations that 
will also be useful in the discussion of reducible representations in the next subsection. These 
are the existence of Cartan (in the semisimple context, maximal Abelian) subalgebras, their con- 
jugacy under the action of the Lie group associated with the algebra, and the fact that any finite- 
dimensional representation, given a choice of Cartan subalgebra (CSA), decomposes into mutually 
orthogonal "weight spaces," which are simultaneously eigenspaces of all CSA elements. The map 
from CSA elements to their eigenvalues on a given weight space is a linear functional on the CSA 
called the "weight" of that weight space. The theorem also uses the observation that the projec- 
tion of the state into the Lie algebra is necessarily a Hermitian element of that algebra, hence 
semisimple (diagonalizable), hence belonging to some CSA, which we call its supporting CSA. 
Frequently, semisimple Lie algebras are presented by giving a Cartan-Weyl basis, consisting of a 
set of commuting, jointly diagonalizable operators that generate a CSA of the algbera, and a set of 
so-called "Weyl operators" that are non-diagonalizable, and act to take a state in one weight space 
to a state in another (or else annihilate it): in physical examples these are often called "raising 
and lowering operators." Normalized states correspond to normalized linear functionals on the 
Lie algebra; when a Cartan-Weyl basis for the algebra is chosen such that the CSA distinguished 
by the basis is the supporting CSA for a given state, the state is zero except on the CSA part of 
the basis. On the CSA, the state is some convex combination of the weights, that is an element 
of the weight poly tope (which is defined as the convex hull of the weights). So it turns out that 
extremal states on the Lie algebra correspond to extremal points of the weight polytope. This 
applies regardless of whether the representation is irreducible or not. For irreducible represen- 
tations (irreps), the extremal points of the weight polytope are also highest-weight states of the 
irrep. Reducible representations are discussed in the next subsection (along with some comments 
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on reductive algebras). 

In preparation for that, we introduce another aspect of standard Lie theory: the Weyl group. 
Besides being able to take any CSA to any other CSA, the Lie group also acts on the weight 
polytope for a given CSA, by reflections in a set of hyperplanes through the origin. The group 
these generate is called the Weyl group. Considered together, the hyperplanes divide the weight 
space into a set of convex cones, sometimes called Weyl chambers, whose points are at the origin, 
and whose union with the hyperplanes is the entire space. Any such cone can be mapped to any 
other via the Weyl group action, and the weight polytope of the representation is the convex hull 
of the Weyl group orbits of the weights in the closure of any single Weyl chamber. 

B. Irreducibly vs reducibly represented Lie algebras 

It is important to realize that the relationships just mentioned between maximal purity, gen- 
eralized coherence, and generalized unentanglement established for a pure state relative to an 
irreducibly represented algebra f) do not automatically extend to the case where f) acts reducibly 
on TL. We will discuss semisimple algebras first and then, because the algebra we use to analyze 
the LMG model is Abelian, the case of reductive algebras. 

If f) is semisimple, a generic finite-dimensional representation of f) may be decomposed as a 
direct sum of irreducible invariant subspaces, H ~ ®{Hi, with each of the Jit being in turn the di- 
rect sum of its weight spaces. Every irrep appearing in the decomposition has a highest (or lowest) 
weight, and for each of these irreps, there is a manifold of GCSs for the irrep constructed as the 
orbit of a highest weight state for that irrep. The weight polytope for the reducible representation 
will be the convex hull of those for all the irreps contained in it. Because of this, the GCSs for these 
irreps will not, in general, all satisfy the extremality property that defines generalized unentangled 
states. This reflects the fact that even for a state belonging to a specific f)-irrep, GE is a property 
which depends in general on how the state relates to the whole representation, not solely the irrep. 
Nor is there necessarily a single weight, for one of the constituent irreps, that generates (as the 
convex hull of the Weyl group orbit) the weight polytope of the reducible representation. Indeed, 
the extremal weights in the weight polytope, which correspond to generalized unentangled states, 
need not all have the same length. Since this squared length is the f)-purity (as defined in Eq. ©) 
of the corresponding state, it is thus no longer the case that all generalized unentangled states have 
maximal Lie-algebraic purity. However, maximal purity remains a sufficient, though no longer a 
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necessary, condition for generalized unentanglement. If the algebra is reductive, the expectations 
of a maximal commutative subalgebra now include ones for the Abelian part of the algebra, i.e. 
operators that commute with the entire algebra. These must be proportional to the identity on each 
irrep, but may have different eigenvalues (possibly degenerate) on different irreps. States on this 
algebra then involve not just weights for the semisimple part of the algebra, but expectation values 
for the Abelian part of the algebra as well. These can distinguish different subsets of the irreps, 
and so irreps whose highest weight for the semisimple part is not extremal for the semisimple part, 
may become extremal (generalized unentangled) in the full reductive algebra. However, maximal 
quadratic purity will remain a sufficient, though in general still not necessary, condition for a state 
being generalized unentangled. 

More intuition about GE, purity, and GCSs may be gained from simple examples. Consider a 
physical system which is composed of two spin-l/2s (namely, two qubits), and let them be labeled 
by A, B, with H = Ha ® Hb — C 4 , and corresponding su(2) generators o^, fff, a G {x, y, z}. 
Consider GE relative to a global representation of su(2), whose total-spin generators are J a = 
a a JrCr a- This representation splits into two irreps, the one-dimensional singlet representation with 
J = and the three-dimensional triplet representation with J = 1. The generalized unentangled 
states relative to this representation of su(2) are those for which there exists an a such that the state 
is a ±1 eigenstate of J a . With respect to the CSA c = { J z }, those are the states | f, t) > 1 1? j)> which 
are also GCSs (with purity equal to 1). No generalized unentangled state is contained in the singlet 
irrep. In particular, neither the spin-zero state in the triplet, nor that which spans the singlet, are 
generalized unentangled (they both have purity equal to 0), nor are they on highest-weight orbits 
(thus GCSs). 

As another example consider a single spin-1 system, whose state space H = C 3 carries an irrep 
of su(2) llln . In this case, for any choice of spin direction (say z) only the J z = ±1 eigenstates 
are generalized coherent. There is also a one-dimensional J z = eigenspace. The maximal-purity 
states are also the highest-weight states; however, the pure J z = eigenstate is not a GCS, has 
zero purity, and is generalized entangled. If, for the same system, a distinguished algebra so (2) 
generated by J z alone is chosen, then the representation reduces as the direct sum of the three 
invariant one-dimensional subspaces corresponding to J z — 1, 0, — 1. In this case, three different 
orbits exist in the representation, each of them consisting of only one state up to phases. However, 
only the states with \J Z \ = 1 are extremal, whereas the state with J z = is not: as one can easily 
verify from the fact that the reduced state is now just the expectation value of J z , an equal mixture 
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of a J z = 1 and a J z = — 1 state has the same reduced state as a J z = state, so the latter remains, 
as in the irreducible case, generalized entangled. 

A generalization of the latter example, which is relevant to the LMG model we will study in 
Section|Vl is the case of a spin- J system with a distinguished Abelian subalgebra generated by J z . 
Again, one can see that only the states with maximal magnitude of J z are generalized unentangled, 
and only they have maximal purity. 

By definition, note that the relative purity and the invariant uncertainty functionals as defined 
in the previous section relate to each other via 

(AF) 2 = (C 2 ) - P h , (8) 

where C 2 denotes the quadratic Casimir invariant of the Lie algebra and P^ is given by Eq. © 
(prior to rescaling). Because, by standard representation theory, C 2 = c e t, with c< £ I within 
each irrep, relative purity and invariant uncertainty essentially provide the same information if f) 
acts irreducibly. This, however, is no longer true in general in the reducible case. In the above 
two-spin- 1/2 example, for instance, the two measures agree on the singlet sector; for triplet states, 
J {J + 1) = 2, thus the invariant uncertainty value is 1 (same as P^) for \J Z \ = 1 (generalized 
unentangled) states, whereas it yields 2 for the (zero-purity) state with J z = in the triplet sector. 

C. Extension to mixed states 

For mixed states on H, the direct generalization of the squared length of the projection onto 
f) as in Eq. © does not give a GE measure with well-defined monotonicity properties under 
appropriate generalizations of the LOCC semigroup of transformations fisl . A proper extension 
of the quadratic purity measure defined in the previous section for pure states to mixed states 
may be naturally obtained via a standard convex roof construction. If p = ^2p s \tp s )(ijj s \, with 

s 

J^Ps — 1 an d Yl pI < 1' m e latter is obtained by calculating the maximum h-purity (minimum 

s s 

entanglement) over all possible convex decompositions {p s , \ip s )} of the density operator p as a 
pure-state ensemble. In general, similarly to what happens for most mixed-state entanglement 
measures, the required extremization makes the resulting quantity very hard to compute. 

While a more expanded discussion of mixed-state GE measures is given in lis! , we focus here 
on applying the notion of GE to characterize QPTs in different lattice systems. Because the latter 
take place in the limit of zero temperature, the ground state of the system may be assumed to be 
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pure under ideal conditions. Accordingly, Eq. © will suffice for our current purposes. 



III. RELATIVE PURITY AS A MEASURE OF ENTANGLEMENT IN DIFFERENT QUANTUM 
SYSTEMS 

We now apply the concept of relative purity to different physical systems in order to understand 
its meaning as a measure of entanglement for pure quantum states. First, we will concentrate 
on spin systems, showing that for particular subsets of observables, the h-purity can be reduced 
to the usual notion of entanglement: the pure quantum states that can be written as a product 
of states of each party will be generalized unentangled. However, for other physically natural 
choices of observable sets, this is no longer the case. Next, we study the h-purity as a measure of 
entanglement for fermionic systems, since this is a good starting point for the analysis of the QPT 
present in the anisotropic XY model in a transverse magnetic field (Section fVTb . In particular, we 
show that if a fermionic state can be represented as a single Slater determinant, it is generalized 
unentangled relative to the Lie algebra u(iV), which is built from bilinear products of fermionic 
operators. These examples illustrate how the concept and measure of GE is applicable to systems 
described by different operator languages, in preparation for the study of QPTs. 

Let us introduce the following representative quantum states for N spins of magnitude S: 



|F£> = \S,S,---,S) , (9) 
i N 

\ w s) = ^]ris,---,s,(s-i) 4 ,s,---,<?>, 



i=l 

2S 



IGHZ^) = - 1 V \S -I, S S -I) , 

where the product state |jSi, S 2 , • ■ ■ , Sn) = \Si)i <8> ® • • • ® \Sn)n, and \Si)i denotes the 
state of the zth party with ^-component of the spin equal to Si (defining the relevant computational 
basis for the ith subsystem). 



A. Two-spin systems 

For simplicity, we begin by studying the GE of a two-qubit system (two spin-l/2s), where the 
most general pure quantum state can be written as \vp) = a||, — |)+c| — \, \)+d\ — ~, — ~), 

with the complex numbers a, b, c, and d satisfying \a\ 2 + \b\ 2 + |c| 2 + \d\ 2 = 1. The traditional 
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measures of pure-state entanglement in this case are well understood, indicating that the Bell states 

| GHZ? ) I37I (and its local spin rotations) are maximally entangled with respect to the local Hilbert 

2 

space decomposition Hi <S>Ti.2- On the other hand, calculating the purity relative to the (irreducible) 
Lie algebra of all local observables f) = su(2) 1 © 5u(2) 2 = {a l a ; i : 1, 2; a — x, y, z} classifies 
the pure two-spin- 1/2 states in the same way as the traditional measures do (see Fig. [T]). Here, the 
operators a\ = a a © 1 and = 11© cr a are the Pauli operators acting on spin 1 and 2, respectively, 
and 

in the basis where | + l/2) 
gives 

where Bell's states are maximally entangled (Pf, = 0) and product states of the form \ip) = 
|0i)i © 1^2)2 (GCSs of the local algebra t) above) are generalized unentangled, with maximum 
purity. Therefore, the normalization factor K = 1/2 may be obtained by setting P^ = 1 in such a 
product state. As explained in Section|nJ P^ is invariant under group operations, i.e., in this case, 
local rotations. Since all GCSs of f} belong to the same orbit generated by the application of group 
operations to a particular product state (a reference state like |~, |) = ||, T))> me Y ai l consistently 
have maximum h -purity (P^ = 1). 

Another important insight may be gained by calculating the purity relative to the algebra of 
all observables for the system, () = su(4) = {o - ^, ct* ® o§; i = 1,2; a, (3 = x,y,z} in this 
case. One finds that any two spin- 1/2 pure state (including Bell's states) is then generalized 
jmentangled (P^ = 1, see also Fig. [T]). This property is a manifestation of the relative nature of 
GE, as considering the set of all observables as being physically accessible is equivalent to not 
making any preferred subsystem decomposition. Accordingly, in this case any pure quantum state 
becomes a GCS of su(4). 

In Fig. [l]we also show the GE for systems of two parties of spin-S* relative to different algebras. 
We observe that the purity reduces again to the traditional concept of entanglement for higher spin 
if it is calculated relative to the (irreducible) Lie algebra of all local observables f) = su(2S' + l)i © 
su(2S' + 1) 2 . For example, if we are interested in distinguishing product states from entangled 
states in a two-spin- Is system, we need to calculate the purity relative to the (irreducible) algebra 



(0 1\ / -i\ I 

, <*x = , o y = , °z = (10) 

= ||) = and |-l/2) = ||) = ^J. In this case, Eq. © simply 

A(I^» = ^E«> 2 ' ( n ) 
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FIG. 1: Purity relative to different possible algebras for a two-spin- 5 system. The quantum states |GHZ|) 
and |F|) are defined in Eqs. @. 

P) = su(3)i © su(3) 2 = {A* <g> I 2 , U 1 <g> A 2 (1 < a < 8)}, where the 3 x 3 Hermitian and traceless 
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matrices Aj are the well known Gell-Mann matrices 112811 : 
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which satisfy Tr^A^] = 5 a ^. In this basis, the computational spin-1 states are represented by 
the 3 -dimensional vectors 



(12) 
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and |-1) = 


















Then, the relative purity for a generic pure state becomes 



a=l i=l 



(13) 

where (A^) denotes the expectation value of A^ in the state In this way, product states like 
|"0) — \4>i)i ® 1^2)2 are generalized unentangled (Pf, = 1) and states like |GHZ^) (and states 
connected through local spin unitary operations), are maximally entangled in this algebra (P, = 0). 

Different results are obtained if the purity is calculated relative to a subalgebra of local observ- 
ables. For example, the two-spin-1 product state |0, 0) = |0) (g) |0) where both spins have zero pro- 
jection along z becomes generalized entangled relative to the (irreducible) algebra su(2)x ©su(2) 2 
of local spin rotations, which is generated by {S l a ; i : 1, 2; a — x, y, z}, the spin-1 angular mo- 
mentum operators S a for each spin being given by 
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Notice that access to local angular momentum observables suffices to operationally characterize 
the system as describable in terms of two spin-1 particles (by imagining, for instance, performing a 
Stern-Gerlach-type of experiment on each particle). Thus, even when a subsystem decomposition 
can be naturally identified from the beginning in this case, states which are manifestly separable 
(unentangled) in the standard sense may exhibit GE (see also Appendix 1X1) . On the other hand, 
this is physically quite natural in the example, since there are no SU(2) x SU(2) group operations 
(local rotations) that are able to transform the state |0, 0) into the unentangled product state |1, 1). 

The examples described in this section together with other examples of states of bipartite quan- 
tum systems are shown in Fig. [T] It is clear that calculating the purity relative to different algebras 
gives information about different types of quantum correlations present in the system. 



B. iV-spin systems 



The traditional concept of pure multipartite entanglement in an N spin-5 1 quantum system 
refers to quantum states that cannot be written as a product of states of each party. The h -purity 
distinguishes pure product states from entangled ones if it is calculated relative to the (irreducible) 

N 

algebra of local observables f} = ® su(2S + l)j (see Appendix lAl. By Eq. ©, the measure P§ is 
invariant under local unitary operations as desired. In particular, the usual concept of entanglement 
in an iV-qubit quantum state (N spin l/2s) can be recovered if the purity is calculated relative to 

N 

the local algebra h = su(2); = [a]., Uy, a\, ■ ■ ■ , erf , a^, a^}, where the Pauli operators a % a 

i=i 

(a = x, y, z) are now 

n factors 



a* a = 1 ® 1 ® • • • <8> a a <g> • ■ ■ ® 1 , (15) 
i th factor 

and the 2x2 matrices a a and 11 are given in Eq. (flOb . Then, the local purity becomes 



a=x,y,z i=l 

where again the normalization factor 1/N is obtained by setting Pf, = 1 in any product state like 

1^) = |0 2 )2 ®- • -<8> \<Pn)n (aGCS in this algebra). With this definition, states like |GHZf ), 

2 

I) — II) ]))/V^] m (with obvious notations), and the cluster states introduced in Ref. 



1 38fl (see also Appendix[B]), will be maximally entangled (Pf, = 0). 

Remarkably, as announced in after some algebraic manipulations (see Appendix Q> one 
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can prove that 

p^m = i - Qim , en) 

where Q is the (pure-state) measure of global entanglement for N spin l/2s systems originally 
introduced by Meyer and Wallach in [39]. A similar relation was independently derived in l^oll . 
See also fl] for additional related considerations. 

N 

In Fig. |2]we display some examples of the purity relative to the local algebra fj = su(2)j for 

-i=i 

a iV spin-5 1 system. We also show the purity relative to the algebra of all observablessu([2S'+l] JV ), 
where any pure quantum state is a GCS, thus generalized unentangled (Pf, = 1). 



C. Purity relative to the u(N) algebra 

We now apply the concept of GE to a physical system consisting of N (spinless) fermion 
modes j, each mode being described in terms of canonical creation and annihilation operators cj, 
Cj respectively, satisfying the following anti-commutation rules: 

{c\, Cj} = 5 itj , {c u Cj} = . (18) 

For instance, different modes could be associated with different sites in a lattice, or to delocalized 
momentum modes related to the spatial modes through a Fourier transform. In general, for any 
N x N unitary matrix U, any transformation Cj t— > J2j ^ij c j ma P s the original modes into another 
possible set of fermionic modes. Using the above commutation relations, one also finds that 

\A c v C k C l\ = S 3k4 c l - S U C k C j ■ ( 19 ) 

Thus, the set of bilinear fermionic operators {etc-,; 1 < j,f < N} provides a realization of 
the unitary Lie algebra u(N) in the 2 N -dimensional Fock space H.Fock of the system. The latter 
is constructed as the direct sum of subspaces H n corresponding to a fixed fermion number n = 
0, . . . , N, with dim(7Y„) = N\/[n\(N — n)\]. For our purposes, it is convenient to express u(N) 
as the linear span of a Hermitian, orthonormal operator basis, which we choose as follows: 



u(iV) 



(c]cj, + c],Cj) with 1 < j < f < N 

i(c] Cj/ - c] lCj ) with 1 < j < f < N , (20) 
V2(c]c 3 - 1/2) with 1 < j < N 

(We use henceforth the notational convention that the large left curly bracket means "is the span 
of"). The action of u(N) on H,Fock is reducible, because any operator in u(N) conserves the 
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total number of fermions n = (X^Li c J c .j)' ^ turns ou t that the irrep decomposition of u(N) 
is identical to the direct sum into fixed-particle-number subspaces Ti n , each irrep thus appearing 
with multiplicity one. 

Using Eq. ©, the h-purity of a generic pure many-fermion state relative to u(N) becomes 

9 N A N 

W)) =mH Wfr + 4s-) 2 - - 4s) 2 l + 77 E(4s- - V2) 2 . (2d 



j<j'=i j=i 

Here, we took K = 2/N, for reasons that will become clear shortly. In this case, the fermionic 
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FIG. 2: Purity relative to different algebras for a iV spin-S" system. The quantum states |GHZg ), | W^), and 
| Fjj) are defined in Eqs. 0. 
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product states (Slater determinants) of the form \<p) = JJ c\\ vac), with |vac) denoting the reference 

i 

state with no fermions and / labelling a particular set of modes, are the GCSs of the u(N) algebra 

nn 

B2L I33H . Because a Slater determinant carries a well defined number of particles, each GCS 
belongs to an irrep space H n for some n, states with different n belonging to different orbits under 
u(N). A fixed GCS has maximum h-purity when compared to any other state within the same 
irrep space. Remarkably, it also turns out that any GCS of h = u(N) gives rise to a reduced 
state which is extremal (thus generalized unentangled) regardless of n, the h-purity assuming the 
same (maximum) value in each irrep. Using this property, the normalization factor K = 2/N was 
calculated by setting Pf, = 1 in an arbitrary Slater determinant. Thus, the purity relative to the 
u(N) algebra is a good measure of entanglement in fermionic systems, in the sense that Pf, = 1 in 
any fermionic product state, and Pf, < 1 for any other state, irrespective of whether the latter has 
a well defined number of fermions or not. Notice that, thanks to the invariance of Pf, under group 
transformations (Eq. ©), the property of a state being generalized unentangled is independent of 
the specific set of modes that is chosen. This is an important difference between our GE and the 
mode entanglement approach 



IV. ENTANGLEMENT AND QUANTUM PHASE TRANSITIONS 

As already mentioned, although many measures of entanglement have been defined in the liter- 
ature, assessing their ability to help us better understand QPTs in quantum systems largely remains 
an open problem. In the following two sections we attempt to characterize the QPTs present in 
the LMG model and in the anisotropic XY model in an external magnetic field through the GE 
notion, relative to a particular subset of observables which will be appropriately chosen in each 
case. Interestingly, for both these models the ground states can be computed exactly by mapping 
the set of observable operators involved in the system Hamiltonian to a new set of operators which 
satisfy the same commutation relations, thus preserving the underlying algebraic structure. In the 
new operator language, the models are seen to contain some symmetries that make them exactly 
solvable, allowing one to obtain the ground state properties in a number of operations that scales 
polynomially with the system size (see also I43I for related discussions). It is possible then to 
understand which quantum correlations give rise to the QPTs in these cases. 

Several issues should be considered when looking for an algebra h of observables that may 
make the corresponding relative purity a good indicator of a QPT. A first relevant observation is 
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that in each of these cases a preferred Lie algebra exists, where the respective ground state would 
have maximum h-purity independently of the interaction strengths in the Hamiltonian. The purity 
relative to such an algebra remains constant, therefore it does not identify the QPT. (In these cases, 
this algebra is in fact the Lie algebra generated by the parametrized family of model Hamiltonians, 
as the parameters are varied.) Thus, one needs to extract a subalgebra relative to which the ground 
state may be generalized entangled, depending on the parameters in the Hamiltonian. A second, 
closely related observation is that the purity must contain information about quantum correlations 
which undergo a qualitative change as the transition point is crossed: thus, the corresponding 
degree of entanglement, as measured by the purity, must depend on the interaction strengths gov- 
erning the phase transition. Finally, whenever a degeneracy of the ground state exists or emerges 
in the thermodynamic limit, a physical requirement is that the purity be the same for all ground 
states. 

Although these restrictions together turn out to be sufficient for choosing the relevant algebra 
of observables in the following two models, they do not provide an unambiguous answer when 
solving a non-integrable model whose exact ground state solution cannot be computed efficiently. 
Typically, in the latter cases the ground states are GCSs of Lie algebras each of whose dimension 
increases exponentially with the system size. Choosing the observable subalgebra that contains 
the proper information on the QPTs (such as information on critical exponents) then becomes, in 
general, a difficult task. 

On the other hand, a concept of generalized mean-field Hamiltonian emerges from these con- 
siderations. Given a Hilbert space H of dimension p N (with p an integer > 1), we will define a 
mean-field Hamiltonian as an operator 



that is an element of an irreducibly represented Lie algebra of Hermitian operators f) = 
{Ax, ■ ■ ■ , Al} whose dimension scales polynomially in N that is, L = vo\y(N). When the ground 
state of such an i/ MF is non-degenerate, it turns out to be a GCS of f) 11511 . while the remaining 
eigenstates (some of which may also be GCSs) and energies can be efficiently computed. The 
connection between Lie-algebraic mean-field Hamiltonians and their efficient solvability deserves 

n 

a careful analysis in its own right, which we will present elsewhere H44I1 . 




(22) 
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V. LIPKIN-MESHKOV-GLICK MODEL 



Originally introduced in the context of nuclear physics 12111 . the Lipkin-Meshov-Glick (LMG) 



model is widely used as a testbed for studying critical phenomena in (pseudo)spin systems W32I1 



This model was shown to be exactly-solvable in I22I1 . In this section, we investigate the critical 
properties of this model by calculating the purity relative to a particular subset of observables, 
which will be chosen by analyzing the classical behavior of the ground state of the system. For 
this purpose, we first need to map the model to a single spin, where it becomes solvable and where 
the standard notion of entanglement is not immediately applicable. 

The model is constructed by considering N fermions distributed in two iV-fold degenerate 
levels (termed upper and lower shells). The latter are separated by an energy gap e, which will 
be set here equal to 1. The quantum number a — ±1 (j or |) labels the level while the quantum 
number k denotes the particular degenerate state in the level (for both shells, k 6 {ki, . . . , k^}). In 
addition, we consider a "monopole-monopole" interaction that scatters pairs of particles between 
the two levels without changing k. The model Hamiltonian may be written as 

H = H + V + W = \Y^ aC la C ka + C lA'a C k'aCka + ^ C lA>a C k'*Cka , (23) 

fc,<r k,k',cr k,k',a 

where a = —a, and the fermionic operators c k(J (c fccr ) create (annihilate) a fermion in the level 
identified by the quantum numbers (k, a) and satisfy the fermionic commutation relations given 
in Section ITlI CI Thus, the interaction V scatters a pair of particles belonging to one of the levels, 
and the interaction W scatters a pair of particles belonging to different levels. Note that the factor 
1 /N must be present in the interaction terms for stability reasons, as the energy per particle must 
be finite in the thermodynamic limit. 

Upon introducing the pseudospin operators 

•4 = 5]4t c *4> ( 24 ) 

k 

J - = Yl c H c k1 > (25) 

k 

J* = g Yl aC ka C ka = g (^T - n l) 7 (26) 

k,a 

which satisfy the su(2) commutation relations of the angular momentum algebra, 

[J z , J±] = ±J ± , (27) 
[J+, J-] = Jz , (28) 
20 



the Hamiltonian of Eq. d23l may be rewritten as 

H = Jz + 2^ (J + + J!) + ^N {J+J - + J - J+) • (29) 
As defined by Eq. (l29l) . if is invariant under the Z 2 inversion symmetry operation K that trans- 
forms (J x , J y , J z ) I— > ( — J x ,—Jy,Jz), and it also commutes with the (Casimir) total angular mo- 
mentum operator J 2 = J 2 + J 2 + J 2 . Therefore, the non-degenerate eigenstates of if are si- 
multaneous eigenstates of both K and J 2 , and they may be obtained by diagonalizing matrices 
of dimension 2 J + 1 (whereby the solubility of the model). Notice that, by definition of J z as in 
Eq. (l26l) . the maximum eigenvalue of J z and J = |J| is N/2. In particular, for a system with 
iV fermions as assumed, both the ground state \g) and first excited state |e) belong to the largest 
possible angular momentum eigenvalue J = N/2 I21I (so-called half-filling configurations); thus, 
they can be computed by diagonalizing a matrix of dimension N + 1. 

The Hamiltonian d29t does not exhibit a QPT for finite N. It is important to remark that some 
critical properties of the LMG model in the thermodynamic limit N — > 00 can be understood by 



using a semiclassical approach 14511 (note that the critical behavior is essentially mean-field): first, 
we replace the angular momentum operators in H/N (with if given in Eq.d29l)) by their classical 
components (Fig. 01); that is 

3 = (J x ,J y ,J z ) — > (J sin 9 cos 0, J sin 9 sin 0, J cos #) , (30) 
if /TV - /i c (j,M) , (3D 

where h c is the resulting classical Hamiltonian and j = J/N, j = 0, . . . , 1/2. In this way, one can 
show that in the thermodynamic limit (see Appendix ITJb 

l im WM. = i im El = min h c (j, 9, 0) , (32) 
jv^oc N at— >oo TV j,e,<f> K ' 

so the ground state energy per particle E g /N can be easily evaluated by minimizing 

h c (j, 9, 0) = j cos 9 + sin 2 9 cos(20) + Wj 2 sin 2 9 . (33) 

As mentioned, the ground and first excited states have maximum angular momentum j = 1/2. 
In Fig. |4]we show the orientation of the angular momentum in the ground states of the classical 
Hamiltonian h c , represented by the vectors J, Ji, and J2, for different values of V and W. When 
A = \V\ — W < 1 we have 9 = it and the classical angular momentum is oriented in the negative 
^-direction. However, when A > 1 we have cos 9 = — A -1 and the classical ground state becomes 
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FIG. 3: Angular momentum coordinates in the three-dimensional space. 

two-fold degenerate (notice that h c is invariant under the transformation <\> i— »■ —(f)). In this region 
and for V < the angular momentum is oriented in the xz plane (0 = 0) while for V > it 
is oriented in the yz plane {<p — ±tt/2). The model has a gauge symmetry in the line V = 0, 
W < — 1, where can take any possible value. 

A. First and second order QPTs, and critical behavior 



Going back to the original Hamiltonian of Eq. d23t . the quantum system undergoes a second 
order QPT at the critical boundary A c = \V C \ — W c = 1, where for A > A c the ground and 
first excited states \g) and |e) become degenerate in the thermodynamic limit and the inversion 
symmetry K breaks. The order parameter is given by the mean number of fermions in the upper 
shell (n-f) = 1/2 + (J z ) /N, which in the thermodynamic limit converges to its classical value, 

1 + cos 6 



lim (n-f) 



(34) 



Obviously, for A < A c we have (n-f) = 0, and (n-f) > otherwise (see Fig|4jl. The critical 
exponents of the order parameter are easily computed by making a Taylor expansion near the 
critical points (A — > 1 + ). Defining the quantities x = V c — V and y = W c — W , we obtain 

(y a - x 13 ) /2 for V > 
(y a + xP)/2 forV < ' 

where the critical exponents are a = 1 and (3 = 1. 

In Fig. 13 we show the exact ground state energy per particle E g /N (with E g = (g\H\g)) as a 
function of V and W in the thermodynamic limit (Eqs. (1321 ). One can see that also in the broken 
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lim (n-t) 



N 2 nd order 




A = |V| - W 

FIG. 4: Representation of the classical ground state of the LMG model. 

symmetry region (A > 1) the system undergoes a first order QPT at V = 0; that is, the first 
derivative of the ground state energy with respect to V is not continuous in this line. 

B. Purity as an indicator of the QPTs in the LMG model 

The standard notion of entanglement is not directly applicable to the LMG model as described 
by Eq. (1291. for this is a single spin system and no physically natural partition into subsystems 
is possible. Therefore, using the h-purity as a measure of entanglement becomes an advantage 
from this point of view, since the latter only depends on a particular subset of observables and no 
partition of the system is necessary. The first required step is the identification of a relevant Lie 
algebra of observables relative to which the purity has to be calculated. 

Since both the ground and first excited states of the quantum LMG model may be understood as 
states of a system carrying total angular momentum J = N/2, a first possible algebra to consider 
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FIG. 5 : Ground state energy per particle in the LMG model. 

is the su(iV + 1) algebra acting on the relevant (N + 1) -dimensional eigenspace. Relative to this 
algebra, \g) is generalized unentangled for arbitrary values of V, W thus the corresponding purity 
remains constant and does not signal the QPTs. However, the family of Hamiltonians d29t do not 
generate this Lie algebra, but rather an 5u(2) algebra, so perhaps su(A r + 1) is not a natural choice 
physically 

Thus a natural choice, suggested by the commutation relationships of Eqs. d27t and d28l) . is 
to study the purity relative to the spin-iV/2 representation of the angular momentum Lie algebra 

f) = SU (2) = {J x , Jy, J Z }\ 



(35) 



where the normalization factor K = N 2 /A is chosen to ensure that the maximum of Pf, is equal 
to 1. With this normalization factor, can be calculated exactly in the thermodynamic limit 
by relying on the semi-classical approach described earlier (see Appendix iDl and Eq. (l30b ^>. For 
V = and arbitrary W > 0, \g) = \J Z = -N/2) which is a GCS of su(2) and has P„ = 1. 
For generic interaction values such that A < 1, the classical angular momentum depicted in Fig. 
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|4]is oriented along the ^-direction and is not degenerate: because (J x ) = (J y ) = 0, only (J z ) 
contributes to Pf,; by recalling that limTv^oo (J^/iV) = —1/2, this gives Pf, = 1, so that as far as 
relative purity is concerned the ground state behaves asymptotically like a coherent state in the 
thermodynamic limit. Physically, this means that with respect to the relevant fluctuations, GCSs of 
su(2) are a good approximation of the quantum ground state for large particle numbers, as is well 
established for this model I^tI . However, in the region A > 1 the ground state (both classical and 
quantum) is two-fold degenerate in the N — »■ oo limit, and the value of Pf, depends in general on 
the particular linear combination of degenerate states. This can be understood from Fig. HJ where 
different linear combinations of the two degenerate vectors Ji and J 2 imply different values of 
(J x ) for V < and different values of (J y ) for V > 0, while (J z ) remains constant. With these 
features, the purity relative to the su(2) algebra will not be a good indicator of the QPT. 

An alternative option is then to look at a subalgebra of su(2). In particular, if we only consider 
the purity relative to the single observable () = so (2) = {J z } (i.e., a particular CSA of su(2)), and 
retain the same normalization as above, we have 

= -^(J*) 2 , (36) 

This new purity will be a good indicator of the QPT, since Pf, = 1 only for A < 1 in the ther- 
modynamic limit, and in addition Pf, does not depend on the particular linear combination of the 
two-fold degenerate states in the region A > 1, where Pf, < 1. Obviously, in this case Pf, is 
straightforwardly related to the order parameter (Eq. (HPT) ): the critical exponents of Pf, — 1 are 
indeed the same (a = 1 and j3 = 1). 

Note that the purity defined by Eq. d36T) does not always take its maximum value for GCSs of 
t) = so (2) (eigenstates of J z ). In the region A < 1 where Pf, = 1, the quantum ground state of the 
LMG model (Eq. (l29t ) does not have a well defined ^-component of angular momentum except 
at V = ([H, J z ] ^ if V ^ 0), thus in general it does not lie on a coherent orbit of this algebra 
for finite N. However, as discussed above, it behaves asymptotically (in the infinite N limit) as a 
GCS (in the sense that Pf, — > 1). Moreover, in Section[n]we showed that for J z -eigenstates with 
eigenvalues | J z \ < N/2, we also obtain Pf, < 1. 

In Fig. |6]we show the behavior of Pf, as a function of the parameters V and W. Interestingly, 
the purity relative to J z is a good indicator not only of the second order QPT but also of the first 
order QPT (the line V = 0, W < -1). 
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FIG. 6: Purity relative to the observable J z in the ground state of the LMG model. 

VI. ANISOTROPIC XY MODEL IN A TRANSVERSE MAGNETIC FIELD 

In this section, we exploit the purity relative to the u(N) algebra (introduced in Section (IIII Q ) 
as a measure able to identify the paramagnetic to ferromagnetic QPT in the anisotropic one- 
dimensional spin- 1/2 XY model in a transverse magnetic field and classify its universality proper- 
ties. 

The model Hamiltonian for a chain of N sites is given by (see Fig|7J 

N N 

H = -a E [C 1 + ^K< +1 + 0- ~ 7K< +1 ] + E < > ( 37 ) 

i=l i=l 

where the operators a' l a (a = x, y, z) are the Pauli spin- 1/2 operators on site i (defined in Eqs. 
(TTOb and (fT5t ). is the parameter one may tune to drive the QPT, and < 7 < 1 is the amount 
of anisotropy in the xy plane. In particular, for 7 = 1 Eq. (l37l) reduces to the Ising model in 
a transverse magnetic field, while for 7 — > the model becomes isotropic. Periodic boundary 
conditions were considered here, that is o^ +iV = a l a , for all i and a. 

When g 3> 1 and 7 = 1 the model is Ising-like. In this limit, the spin-spin interactions are the 
dominant contribution to the Hamiltonian d37t . and the ground state becomes degenerate in the 
thermodynamic limit, exhibiting ferromagnetic long-range order correlations in the x direction: 
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FIG. 7: Anisotropic one-dimensional XY model in an external transverse magnetic field B. 



M% = limAr_ >00 (o"^,(T^ 2 ) > 0, where M x is the magnetization in the x-direction. In the opposite 
limit where g — > 0, the external magnetic field becomes important, the spins tend to align in the z 
direction, and the magnetization in the x direction vanishes: = Ym\ N ^ 00 (a].(Jx^ 2 ) = 0. Thus, 



in the thermodynamic limit the model is subject to a paramagnetic-to-ferromagnetic second order 
QPT at a critical point g c that will be determined later, with critical behavior belonging to the 2-D 
Ising model universality class. 



This model can be exactly solved using the Jordan-Wigner transformation [48], which maps 
the Pauli (spin 1/2) algebra into the canonical fermion algebra through 



j'-i 

n< 



+ ' 



(38) 



where the fermionic operators cj (cj) have been introduced in Section HlTCl and o\ = (cr£ + ia J y )/2 
is the raising spin operator. 

In order to find the exact ground state, we first need to write the Hamiltonian given in Eq. (l37l 
in terms of these fermionic operators, 



N-l 



H = -2gJ2(4c l+ i + 744+1 + h.c.) + 2gK(J N c 1 + ^ N c\ + h.c.) + 2N 



(39) 



i=l 



N N 

where K = J| (— a J z ) is an operator that commutes with the Hamiltonian, and N = ^ cjcj is the 

j=i i=i 
total number operator (here, we choose iV to be even). Then, the eigenvalue of K is a good quan- 



tum number, and noticing that K 



JnN 



we obtain K 



■1) whenever the (non-degenerate) 
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eigenstate of if is a linear combination of states with an even (odd) number of fermions. In partic- 
ular, the numerical solution of this model in finite systems (with N even) indicates that the ground 
state has eigenvalue K = + 1, implying anti-periodic boundary conditions in Eq. d3*9"l ). 

The second step is to re-write the Hamiltonian in terms of the fermionic operators c\ (c k ), 
defined by the Fourier transform of the operators cj (c ): 

N 



TtX^c), (40) 
i=i 

where the set V of possible k is determined by the anti-periodic boundary conditions in the 
fermionic operators: V = V + + = [±^, ±7^- , • ■ ■ , db ^^ 1 ^ ]. Therefore, we rewrite the Hamil- 
tonian as 

H + N = -2J2(~1 + 29 cos k)& k + i gi sin k(dl k c{ + c_ k c k ) . (41) 

keV 

The third and final step is to diagonalize the Hamiltonian of Eq. (BIT) using the Bogoliubov 
canonical transformation 

Ik = u k c k - iv k c ] _ k 
ll k = u k cl k - iv k c k 
where the real coefficients u k and v k satisfy the relations 

u k = u_ k , v k = -v_ k , and u 2 k + v 2 k = l , (42) 



where 

?/,,„ = ms I 

2 ) ' V 2 



u k = cos ( ^) , v k = sin (-^ ) , (43) 



with 4> k given by 



. , . 207 sin k ^ 

tan fc = ~~TZo 1 ■ (44) 

— 1 + 2(7 cos k 

In this way, the quasiparticle creation and annihilation operators 7^ and 7 fc , satisfy the canonical 
fermionic anti-commutation relations of Eq. (fTSt . and the Hamiltonian may be finally rewritten as 

H = J2Ulhk~im, (45) 

keV 

where C, k = 2a/ (—1 + 2(7 cos k) 2 + 4g 2/ y 2 sin 2 k is the quasiparticle energy. Since in general 
£ k > 0, the ground state is the quantum state with no quasiparticles (BCS state 14911 ). such that 
7fc|BCS) = 0. Thus, one finds 

I BCS) = J] («* + ™ k clcl k ) |vac) , (46) 

keV+ 
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where |vac) is the state with no fermions (cfc|vac) = 0). 

Excited states with an even number of fermions (K = +1) can be obtained applying pairs 
of quasiparticle creation operators j k to the |BCS) state. However, one should be more rig- 
orous when obtaining excited states with an odd number of particles, since K — — 1 implies 
periodic boundary conditions in Eq. (1391 . and the new set of possible k's (wave vectors) is 
V = [-7T, 0, ^, ... , ^y^] (different of V). 

A. QPT and critical point 

In Fig. [8] we show the order parameter M x = \imN-^oo(o'l<Jx^ 2 ) as a function of g in the 
thermodynamic limit and for different anisotropics 7 |4]. We observe that = for g < g c and 
M x 7^ for g > g c , so the critical point is located at g c = 1/2, regardless of the value of 7. The 
value of g c can also be obtained by setting = in Eq. d45t . where the gap vanishes. 

Notice that the Hamiltonian of Eq. (l37l) is invariant under the transformation that maps 
(a x ; a 3 y \ a k z ) t— > (— a x ; — a^; a z ) (Z 2 symmetry), implying that (cr x ) = for all g. However, since in 
the thermodynamic limit the ground state becomes two-fold degenerate, for g > g c , it is possible 
to build up a ground state where the discrete Z 2 symmetry is broken, i.e. (<j x ) 7^ 0. This statement 
can be easily understood if we consider the case of 7 = 1, where for < g < g c the ground state 
has no magnetization in the x direction: For g = 0, the spins align with the magnetic field, while 
an infinitesimal spin interaction disorders the system and M x = 0. On the other hand, for g — > 00 
the states \g x ) = ^[|-v ■,->) + |<-, ■■•,<-)] and I02) = 75 [H, ' • ' , ^> - \*~, ■ • ' , <-)], with 
|— = ^[|T) + ID] an d l^) = 75[|T) — II)] become degenerate in the thermodynamic limit, and 
a ground state with (a l x ) 7^ can be constructed from a linear combination. 

Remarkably, this paramagnetic-to-ferromagnetic QPT does not exist in the isotropic limit 
(7 = 0). In this case, the Hamiltonian of Eq. dT7t has a continuous u(l) symmetry; that is, it 
is invariant under any z rotation of the form exp[i9 J2j a i\- Since the model is one-dimensional, 
this symmetry cannot be spontaneously broken, regardless of the magnitude of the coupling con- 
stants. Nevertheless, a simple calculation of the ground state energy indicates a divergence in its 
second derivative at the critical point g c = 1/2, thus, a second order non-broken symmetry QPT. 
For g < g c all the spins (in the ground state) are aligned with the external magnetic field, with 
total magnetization in the z direction M z = Yl,j{ a i) = ~ ^> an ^ the quantum phase is gapped. 
For g > g c , the total magnetization in the z direction is M z > —N, the gap vanishes, and the 
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FIG. 8: Order parameter M% in the thermodynamic limit as a function of g for different anisotropics 7. The 
critical point is at g c = 1/2. 

quantum phase becomes critical (i.e., the spin-spin correlation functions decay with a power law), 
with an emergent u(l) gauge symmetry 111 711 . Then, in terms of fermionic operators (Eq. (1391)). 
an insulator-metal (or superfiuid) like second order QPT exists at g c for the isotropic case, with no 
symmetry order parameter. It is a Lifshitz transition. 

B. u(iV)-purity in the BCS state, and critical behavior 

The I BCS) state of Eq. (|4"6T ) is a GCS of the algebra of observables fj = so(2N), spanned by an 
orthonormal Hermitian basis which is constructed by adjoining to the basis of u(iV) given in Eq. 
(l20l the following set r of number- non-conserving fermionic operators: 



Then, the | BCS) state is generalized unentangled with respect to the so(2iV) algebra and its purity 
Pf) (Eq. ©) contains no information about the phase transition: Pf, = 1 V<7, 7. Therefore, in 





(47) 
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order to characterize the QPT we need to look at the possible subalgebras of so(2iV). A natural 
choice is to restrict to operators which preserve the total fermion number that is, to consider the 
u(N) algebra defined in Section UTTC! relative to which the |BCS) state may become generalized 
entangled. (Note that as mentioned in Section lTlICl the u(N) algebra can also be written in terms 
of the fermionic operators ct and c fe , with k belonging to the set V .) 

In the | BCS) state, (c\c k ,) ^ only if k = k' , thus using Eq. (|2"TT) the purity relative to 
() = u(N) is: 

A(|BCS» = 1 £<cfo - 1/2) 2 = ± ~ V2) 2 , (48) 

keV keV 

where the coefficients Vk can be obtained from Eqs. (l43l) and (l44b . In particular, for g = the 
spins are aligned with the magnetic field and the fully polarized | BCS) g=0 = |i, !,•••, 1) state is 
generalized unentangled in this limit (a GCS of u(N) with p, = 1). In the thermodynamic limit, 
the purity relative to the u(N) algebra can be obtained by integrating Eq. d48l) : 

P„(|BCS» = ^ j(v 2 k -l/2) 2 dk, (49) 
o 

leading to the following result: 

ifg< 1/2 



ndBcs)) 



(50) 



Vl-4 5 2(l-7 2 ). 

TT^ if y > 1/2 

Although this function is continuous, its derivative is not and has a drastic change at g = 1/2, 
where the QPT occurs. Moreover, Pf, is minimum for g > 1/2 implying maximum entanglement 
at the transition point and in the ordered (ferromagnetic) phase. Remarkably, for g > 1/2 and 
N — > oo, where the ground state of the anisotropic XY model in a transverse magnetic field is 
two-fold degenerate, Pf, remains invariant for arbitrary linear combinations of the two degenerate 
states. 

As defined, for large g the purity Pf, approaches a constant value which depends on 7. It 
is convenient to remove such dependence in the ordered phase by introducing a new quantity 
= Pf, — (shifted purity). We thus obtain 

ifg< 1/2 

(51) 



n(l BCS » = 



v/l-V(l- 7 2) 

if 5 > 1/2 



The new function P^ behaves like a disorder parameter for the system, being zero in the ferro- 
magnetic (ordered) phase and different from zero in the paramagnetic (ordered) one. The behavior 
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of as a function of g in the thermodynamic limit is depicted in Fig. |9]for different values of 7. 
In the special case of the Ising model in a transverse magnetic field (7 = 1), one has the simple 
behavior P^ = 1/2 - 2g 2 for g < 1/2 and P^ = if g > 1/2. 
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FIG. 9: Shifted purity P^f N \ of the |BCS) as a function of g for different anisotropics 7, Eq. d5TT i. P^ N ^ 
behaves like a disorder parameter for this model, sharply identifying the QPT at g c = 1/2. 

The critical behavior of the system is characterized by a power-law divergence of the correla- 
tion length e, which is defined such that for g < 1/2, lim|i_ 3 -|_ +00 |(<t*ct^)| ~ exp(— Ji^l). Thus, 
e — > 00 signals the emergence of long-range correlations in the ordered region g > 1/2. Near 
the critical point (g — > 1/2 - ) the correlation length behaves as e ~ (g c — g)~ v , where v is a crit- 
ical exponent and the value v — 1 corresponds to the Ising universality class. Let the parameter 
A 2 = e _1 / £ . The fact that the purity contains information about the critical properties of the model 
follows from the possibility of expressing P^ for g < 1/2 as a function of the correlation length, 



i^(|BCS» = 



7 



(52) 



2<?A 2 (l- 7 )-l 

where a known relation between g, 7, and A 2 has been exploited fl]. Performing a Taylor expan- 
sion of Eq. (I52"l) in the region <? — »■ l/2~, we obtain P^ ~ 2(g c — g) u fj with v — 1 and 7 > (Fig 
TOl) . Thus, the name disorder parameter for P^ is consistent. 
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FIG. 10: Scaling properties of the disorder parameter for anisotropy 7 = 1. The exponent v = 1 belongs to 
the Ising universality class. 

Some physical insight in the meaning of the ground-state purity may be gained by noting that 
Eq. (l48t can be written in terms of the fluctuations of the total fermion operator N 

P b (|BCS)) = l-|((iV 2 )-(iV) 2 ) . (53) 

where the |BCS) -property (c\c k >) = Sk,k'vl has been used. In general, the purity relative to a 
given algebra can be written in terms of fluctuations of observables II 1 311 . Since fluctuations of 
observables are at the root of QPTs it is not surprising that this quantity succeeds at identifying 
the critical point. Interestingly, by recalling that P S0 (2tv)(|BCS)) = 1, the u(N) -purity can also be 
formally expressed as 

P uW (|BCS)) = l-^(A Q ) 2 , (54) 

A a ex 

where the sum only extends to the non-number-conserving so(2iV)-generators belonging to the 
set r specified in Eq. (l47b. Thus, the purity is entirely contributed by expectations of operators 
connecting different u(A r )-irreps, the net effect of correlating representations with a different par- 
ticle number resulting in the fluctuation of a single operator, given by N = J2k m ^8- EH 
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we show the probability fi(n) of having n fermions in a chain of N = 400 sites for 7 
observe that for g > 1/2 the fluctuations remain almost constant, and so does the purity. 
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FIG. 1 1 : Distribution of the fermion number in the | BCS) state for a chain of N = 400 sites and anisotropy 
7 = 1. 

Again, the isotropic case (7 = 0) is particular in the sense that P^ = 1 (or P^ = 0, see Fig. 
|9l), without identifying the corresponding metal-insulator QPT. The reason is that in this limit, the 
Hamiltonian of Eq. d39t contains only fermionic operators that preserve the number of particles 
(i.e., H G u(N)), and the ground state of the system is always a GCS of the u(N) algebra. 
Therefore, in order to obtain information about this QPT, one should look into algebras other than 
u(iV), relative to which the ground state is generalized entangled. For example, in Sec. IVIDI 
we study the purity relative to the local algebra of observables and in Fig. [13] we show that it 
succesfully identifies the QPT in the isotropic case, being maximum for g < g c (thus implying 
generalized unentanglement). 

C. Comparison with concurrence 

As mentioned, the critical behavior of the XY model in a transverse field has also been investi- 
gated by looking at various quantities related to the concurrence, which is intrinsically a measure 
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of bipartite entanglement. For a generic mixed state p of two qubits, the latter is calculated as J14I1 

C(p) = max{Ai - A 2 - A 3 - A 4 , 0} , 

where Ai > . . . > A4 are the square roots of the eigenvalues of the matrix R = pp and p = a y ® 
<y y p*<y y <8> cr y . The concurrence for the reduced density operator p^ m of two nearest-neighbor qubits 
(\£ — m\ = 1) and next-nearest-neighbor (\£ — m\ = 2) qubits on a lattice has been investigated 
in detail in Ref. [8]. Since, thanks to translational invariance, p^ m depends on the qubit indexes 
only via their distance, we will use the notation C(l), C(2) for the resulting quantities as in Jsl. 
While the results reported in the above work nicely agree with the scaling behavior expected for 
this model, the emerging picture based on concurrence cannot be regarded as fully satisfactory. As 
stressed in J^], the entanglement as quantified by the nearest-neighbor concurrence is not directly 
an indicator of the QPT in this model, showing maximum entanglement at a point which is not 
related to the QPT. However, the derivative dC(l) / dg of the concurrence with respect to the spin- 
spin coupling parameter can be seen to diverge logarithmically at the critical point for 7 > 0, and 
with a power law for the isotropic case [Fig. (fT2b l. identifying the critical point in this model. 
Such a divergence is not found when analyzing, at the isotropic point, other QPTs in models of 
interest, like the one-dimensional anisotropic Heisenberg chain (see, for instance, [51]). Therefore, 
it suggests that the identification of a critical point using concurrence could be a hard task in 
general. 



D. Purity of the BCS state relative to the local algebra 



Finally, we have also investigated the behavior of the purity of the | BCS) state relative to the 

N 

algebra of local observables f) = ® su(2);. Using Eq. (fToT) . this is physically related to the total 

i=X 

magnetization M 2 Z along z. The resulting behavior is plotted in Fig. [T^Jas a function of g and 
7. As explained in Section IIII Bl this is a measure of the usual notion of entanglement in the 
iV-spin-l/2s system. In particular, the | BCS) state is unentangled for g — > (where | BCS) ~ 
|— §)i® • • -<8> |— |) at), thus Pj, — >• 1 in this limit. Moreover, for g — > 00 we have | BCS) ~ (GHZ^) 
(up to local rotations), thus | BCS) becomes maximally entangled, and — > 0. 

N 

Compared to the purity relative to the u(N) algebra, the purity relative to h = ® su(2)j is 
not as good an indicator of the phase transition when 7 > 0, in the sense that it does not present 
any drastic change in its behavior. However, its derivative with respect to the spin- spin coupling 
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FIG. 12: Nearest-neighbor concurrence and its derivative for the |BCS) state as a function of g in the 
isotropic XY model, 7 = 0. Both curves correspond to the exact solution in the thermodynamic limit. The 
value of dC(l)/dg below g c is also zero as C(l) (not shown). 

parameter diverges at the critical point in this model [Fig. \T%&. Only in the isotropic case (7 = 0) 
the purity relative to the local algebra presents a drastic change at the critical point (see Fig. [T3l . 
In this case, the operator M z = a i) f° r 9 ~~ * dt scales as 

M z + l~(g- g c y (55) 

with the exponent being x = 1/2- On the other hand, this exponent can also be obtained from the 
purity relative to the local algebra, in the same limit: 

l-I\~{g- g c ) x . (56) 
Therefore, this measure of entanglement is also a good indicator of the QPT for the isotropic case. 

VII. CONCLUSIONS 

In this paper, we have explored the usefulness of generalized entanglement (GE) for character- 
izing the broken (and one example of non-broken) symmetry quantum phase transitions (QPTs) 
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JV 

FIG. 13: Purity of the |BCS) state relative to the local algebra su(2)j, as a function of g for different 

i=i 

anisotropies 7 (g c = 1/2). The number of sites N = 400 as in Fig.fTTl 

present in different lattice systems. As we focused on situations where the physically relevant 
observables form a Lie algebra, a natural GE measure provided by the relative purity of a state 
relative to the algebra has been used to identify and characterize these transitions. 

In Sections !!!! Al and flll BL we showed using several illustrative examples how the concept of Im- 
purity can be useful for different spin systems, by encompassing the usual notion of entanglement 
if the family of all local observables is distinguished. In addition, the possibility to directly apply 
the GE notion to arbitrary quantum systems, including indistinguishable particles, was explicitly 
shown in Section ITlI CI using fermionic systems as a relevant case study. Depending on the sub- 
set of observables chosen, the h-purity contains information about different n-body correlations 
present in the quantum state, allowing for a more general and complete characterization of en- 
tanglement. Finally, in Sections |V] and IVj we showed that the h-purity successfully distinguishes 
between the different phases present in two lattice systems, where the critical points are charac- 
terized by a broken symmetry (or non-broken symmetry in the case of the isotropic XY model 
in an external magnetic field) and the usual notion of entanglement cannot be straightforwardly 
applied. As also discussed in Section |lVj the most critical step is to determine which subset of 
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FIG. 14: Derivative of the purity of the | BCS) state relative to the local algebra as a function of g for 7 = 0.5 
and different lattice sizes. 

observables may be relevant in each case, since the h-purity must contain information about the 
quantum correlations that play a dominant role in the QPT. In particular, the ground state of the 
two models we considered can be exactly calculated and the relevant quantum correlations in the 
different phases are well understood, thus choosing this subset of observables becomes relatively 
easy. 

Applying these concepts to a more general case, where the ground state of the system cannot be 
exactly computed, can be done in principle by following the same strategy. However, determining 
in a systematic way the minimal subset of observables fj whose purity is able to signal and char- 
acterize the QPT, thereby providing the relevant correlations, requires an elaborate analysis. Even 
more interesting, perhaps, is the open question of finding the minimal number of GE measures, 
possibly including measures of GE relative to different observable sets, needed to unambiguously 
characterize the universality class of a transition, obtaining all of its critical exponents. Finally, a 
fascinating direction for further investigation is to explore the significance of the GE notion within 
topological quantum-information settings I52I and to understand what generalizations might be 
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needed to handle topological QPTs. 
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APPENDIX A: SEPARABILITY, GENERALIZED UNENTANGLEMENT, AND LOCAL PURI- 
TIES 

Given a quantum system S whose states belong to a Hilbert space H of dimension 
dim(7i) = d, the purity relative to the (real) Lie algebra of all traceless observables f) = su(d) 
spanned by an orthogonal, commonly normalized Hermitian basis {Ai ■ ■ ■ A L }, L = d 2 — 1, is, 
according to Eq. ©, given by: 

L 

A(hA» = K]T(A Q > 2 . (Ai) 

a=l 

The normalization factor K depends on d and is determined so that the maximum purity value is 1 . 
If Jx{A a Aj3) = Sa t/ 3 (as for the standard spin-1 Gell-Mann matrices), then K = d/(d — 1), whereas 
in the case Tx(A a Ap) = d5 a ^ (as for ordinary spin-1/2 Pauli matrices), K — l/(d — 1). Recall 
that any quantum state |V>) G H can be obtained by applying a group operator U to a reference 
state | ref) (a highest or lowest weight state of su(gQ); that is 

|^> = [/|ref>, (A2) 

with U = e t ^' c ' taAa , and t a real numbers. Therefore, any quantum state is a GCS of su(d), 
thus generalized unentangled relative to the algebra of all observables: P^tp)) = 1 for all 

Let now assume that S is composed of N distinguishable susbsytems, corresponding to a fac- 
torization H = <S>f =1 'Hp with dim(7Yj) = d p d = Ylj dj. Then the set of all local observables 
on S becomes f) = f)/ oc = ® . su(gL,). An orthonormal basis which is suitable for calculating the 
local purity Pf, may be obtained by considering a collection of orthonormal bases {A J ai ■ ■ ■ A J aL }, 
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Lj — <B — 1, each acting on the jth subsystem that is, 



n factors 



Ai 3 = I 1 ® ll 2 ® • • ■ ® A a . ® • • • ® 1* , (A3) 

j th factor 

where F = 1/ y^. Then for any pure state \vp) G 7i one may write 

N Lj 

(A4) 



w»=K # x;[E<4 f > s 



j=l aj =l 

By letting hj = span{A Qj .} be the Lie algebra of traceless Hermitian operators acting on Hj alone, 
the above equation also is naturally rewritten as 

N 



n(i^)) = K / E^n J (i^)), K^-f-, ( A5) 

j=l ■? 3 

The h., -purity P^. may be simply related to the conventional subsystem purity. Let pj = 
Tr i7 i J ({|^)('0|}) be the reduced density operator describing the state of the jth subsystem. Be- 
cause the latter can be represented as 

Pi = J+Y< ( A «M«i = £ (Ai.)A aj , (A6) 

3 a,=l ct,=l 



one can also equivalently express Eq. (IA4I) as 

W» = K'X;[Tr^-i 



(A7) 



that is, P^.(|^)) = (djTrp| — l)/(dj — 1). Clearly, the maximum value of either Eqs. (IA5I) or (IA7I) 
will be attained when, and only when, each of the conventional purities TrpJ = 1 <-> = 1 for 
all j, which allows determining the K' -normalization factor as 

K ' = v^"r = r, ^ i = ' \ ■ < A8 > 

Ei4 /v(i-iE4) 

Accordingly, 

P„ !oo (|^» = max = 1 <- |^> = |0i) ® • • • ® |<M , (A9) 

and the equivalence with the standard notions of separability and entanglement are recovered. 
Note that for the case of iV qubits considered in Section IIIB, the above value simplifies to K' = 
2/N which in turn gives the purity expression of Eq. (fToT) once the standard unnormalized Pauli 
matrices are used (A 3 a . = ct£ / y/2, thus removing the overall factor 2). 
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APPENDIX B: CLUSTER AND VALENCE BOND SOLID STATES ARE MAXIMALLY EN- 
TANGLED 



In Ref. 113 811 . Briegel and Raussendorf introduced the so-called cluster states for a system of N 



qubits in D space dimensions which, in the computational basis, are expressed as 

\®)c = 2*72 ® ( lt>i II ^' +7) + li>; ) > ( B1 ) 

jec \ 7 er / 
where C defines the cluster (C C Z D ) and 7 denotes some nearest neighbor qubits in the cluster: 
r = {1} for D = 1, T = {(1, 0), (0, 1)} for D = 2, T = {(1, 0, 0), (0, 1, 0), (0, 0, 1)} for D = 3, 
etc. We consider ai J+7 ' ) = 1 when j ' + 7 is not in C. 

The usual notion of entanglement, as applied to a cluster state, is recovered when the f)-purity 

is calculated relative to the local algebra h = su(2)j (see Appendix lAl. For this purpose, we 

j&c 

first calculate the expectation values (cr J a )c, with a = x,y, z. One can immediately realize that 
{a y )c = 0, Vj, since a 3 y is an Hermitian operator (i.e., (a j y ) G R) that acting on the j-th qubit's 
state (in the natural basis) introduces a phase factor ±i, and the coefficients of Eq. (IB II) are all 
real. Moreover, (er^)c = 0, Vj, since the weight of every state of the natural basis is the same in 
Eq. (IB It . In other words, we have a linear combination of basis states where each single qubit 
has the same probability of pointing up or down. Finally, one can also prove that {o{) c = 0, Vj. 
This can be done by using the eigenvalue equations Kj\$)c = ±|$)c for the family of operators 
Kj = a i ll 7 ^ where f = T [j —T denotes the set of all nearest-neighbor qubits to the j-th 

qubit. Therefore, {a{) c = ±(<riK s ) c = MU Again, since Eq. (IB 1 1> is a combination 

7er 

of all the states of the computational basis with the same probability, we obtain (a J x )c = 0. In this 
way, the h-purity (Eq. dT6b l is Pf, = 0, and the cluster states are maximally entangled relative to 
the local set h = su(2)j. 

Another important class of spin states is the one defined by the so-called Valence Bond Solid 
(VBS) states. These states have been introduced in the context of Heisenberg-like magnets, and 
have been recently revisited in the context of quantum computation H50I1 . Their general form is 

i$)vbs = n (°h - & i a i) M < b2) 

where represent nearest-neighbor bonds of a D-dimensional lattice of coordination z, aj 
and bj are Schwinger-Wigner boson (creation) operators on site j whose relation to su(2) spin-^ 
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generators is 

Si = \{a% + b] aj ), SI = !(<,}&, - bfa), Si = \{a)a 3 - (B3) 

with the constraint aja^- + bjbj = 2S, and M = 2S/z. M being an integer makes the possible 
values of S to depend upon the connectivity of the lattice, which is defined by z. 

We start by showing that the bond operators a\ b^ — b\aj are invariant under global spin rotations. 
The Schwinger-Wigner boson operators transform as vectors for su(2) rotations 




cos §2. gi(0 3 +0i)/2 gj n §2. e i(03-0l)/2 

- sin ^ e -*(03-0i)/2 cos §2 e -i(e 3 +e 1 )/2 




(B4) 



under an arbitrary spin rotation on lattice site j, defined by 

U j = e WlSi e id2S y e^ si , UjU} = U]U j = t . (B5) 
Then, we can use this result to prove that 

U^iajb] - b\a\)U\u) = a\b] - b\a) , (B6) 

implying, for W = ^ llj, 

U j |$) VBS = |$>VBS • (B7) 

Therefore, |$)vbs belongs to the singlet irrep of the total spin J a = J2j $1 (i- e -' (^)vbs = 0). 

We want to show now that (S^)vbs = 0, Vj. We first observe that (S'f)vBS = 0, Vj, because 
the transformation that maps a)- \— > bj and bj i— > — at (i.e., a global spin rotation about the y-axis, 
setting 6i = 6> 3 = and 9 2 — tt in Eq. (IB4I) ) implies (aja^vBS = (^j^)vBS- Then, from the 
invariance under global spin rotations and the singlet nature of |$)vbs> we obtain (S 3 x )y B s = = 

(S^)vbs, Vj. In other words, the purity relative to the algebra f) = 0su(2)j vanishes, meaning 

j 

that |$)vbs is maximally generalized entangled relative to this algebra. 

However, in order to make contact with the standard notion of entanglement (Appendix |A]) we 

need to address the GE relative to the algebra f) = 0su(2S' + that is, relative to the set of 

j 

all local observables. For simplicity, we only discuss the ID case for S = 1 (i.e., M = 1 in Eq. 
SB21 ) but the reader could use the same techniques to obtain results in higher D dimensions and 
spin magnitude S. 

The algebra f> = 0su(3),- = {Sp, 

j 

K*A\ = *0 (v» - v (B8) 
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can be written in terms of the su(2) generators as 







°oo — 


^ f en 2 cj ^2 + 1) 
3 {bz) ' 611 " 2 


1 

~ 3 




^[5i+{4,^}] 




°01 — 


J^[£+{si,^}] 




°20 — 


-L[sL-{sL,s{}] 

— v ^ 




°02 — 






°12 — 




-1, 


°21 — 


1 2 1 • 2 

— {iS 1 ^, S^} + (•S'x) 2 + 2 (^) 2 


-1, 



(B9) 



with S*4 = ± iS^. From the spin-rotational invariance of the state |$)vbs we g et )vbs — 

((^) 2 )vbs = ((^') 2 )vbs = and > sin ce 5 = 1, we obtain (^ ) VBS = (^i)vbs = 0. 

Moreover, the spin-rotational invariance also implies that (S^S 3 a ,) vbs remains the same constant 
Va ^ a'. Then, for example, applying a global n -rotation about the y-axis to |$)vbs (i-e., the 
operation that maps S J Z 1— > — and h-> Sp we obtain (S^S J z )y B s — —(S^S J z )y B s — 0, hence, 
(«S^}vbs = 0. Therefore, the state |$)vbs (5 1 = 1, M = 1) is maximally entangled when using the 
standard notion of entanglement {P^ = 0, for the algebra of all local observables h = su(3) 3 -). 

APPENDIX C: RELATION BETWEEN PURITY IN THE LOCAL ALGEBRA AND THE 
MEYER- WALLACH MEASURE OF ENTANGLEMENT 



In Ref. 113 911 . Meyer and Wallach define a measure of entanglement Q for pure states of qubit 
systems, that is invariant under local unitary operations (local rotations). For this purpose, they 
first define the mapping lj(b) acting on product states as 

lj{b)\bi, ■ ■ - ,b N ) = Sbb^bi, - ■ -,bj, ■ ■ -,b N ), (CI) 

where b and bj are either the states ||) or |— |), and bj denotes the absence of the j-th qubit. On 
the other hand, any iV-qubits pure quantum state can be written in the natural basis (^-component 
of the spin equal to ±|) as 

2JV-1 

M = E [d\\)i + K\-\)]\<t>i), (C2) 



i=l 
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where gj and h\ are complex coefficients, and the orthonormal states \(j>i) of N — 1 qubits (absence 
of the j-th qubit) are also written in the natural basis. Therefore, the action of lj(b) on is 



2-iV-i 

8=1 
2 iV-l 

i=l 



Then, they define the entanglement Q(\ip)) as 



4 N 1 1 

^l^)) = ^E D ( / ^2 )l?/;) ' /j( ~2 )l?/;) 



(C3) 



where the distance between two quantum states \u) = ^2 an d \v) = ^2 v i\4>i) 1S 

D(u, v) = - E \ u i v i ~ 



UjVi\ 



(C4) 



i -j 



Therefore, 



(C5) 



where * denotes complex conjugate. After some simple calculations we obtain the following 
relations 



E JM) 



i=i 



2 

1 - ^ 



(C6) 
(C7) 
(C8) 



i=l 



and the distance becomes D(/j(i)|^), ^-(-|) |^)) = J[l - (^) 2 - (a£) 2 - (a^) 2 ]. Since Q(|^)) 
contains a sum over all qubits (see Eq. (IC3I) ). we finally obtain 



N 



Q(m) = i - ^ E I <^> 2 + + K> 2 = 1 - p *(i^» 



(C9) 



N 



where is the purity relative to the local algebra h; oc = su(2)j defined in Section lTlIBI 



j'=i 
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APPENDIX D: CLASSICAL LIMIT IN THE LMG MODEL 



As we mentioned in Section|Vl some critical properties of the LMG, such as the order parameter 
or the ground state energy per particle in the thermodynamic limit, may be obtained using a semi- 



classical approach. In this section we sketc 



i a rough analysis of why such approximation is valid 



(for a more extensive analysis, see Ref. 114511 ') . 
We first define the collective operators 

N 

E {°,0')=^ C k0 C k<T<l ( D1 ) 
fc=l 

where a, a' =f or j and the fermionic operators c ka . (c ka ) have been defined in Section |Vj The 
collective operators satisfy the u(2) commutation relations (Section Till CI ): that is 

[E{a,a>), -EVV")] — $ct'<t"E(ct,<t'") — 5 acT '» E( a 'r >a 'y (D2) 

If the number of degenerate levels N is very large, it is useful to define the intensive collective 
operators Eu a <) = E/ a>a r\/N, with commutation relations 

E( a ,a>), E( a iiyii) 



y8a'a"E( aia iii) — Saa'i'E^ir ^j . (D3) 

Therefore, the intensive collective operators commute in the limit N — > oo, they are effectively 
classical and can be simultaneously diagonalized. Similarly, the intensive angular momentum 
operators J x /N — (E^^+E^))/2, J y /N — (E^^ — E^)/2i, andJ z /N — (E^^) — E^ t i))/2 
(with J a defined in Eqs. (l24b . (1251) . and d26j)) commute with each other in the thermodynamic limit, 
so they can be thought of as the angular momentum operators of a classical system. 

Since the intensive LMG Hamiltonian H/N, with H given in Eq. (I29L can be written in terms 
of the intensive angular momentum operators, it can be regarded as the Hamiltonian describing a 
classical system. The ground state of the LMG model \g) is then an eigenstate of such intensive 
operators when N — > oo: (J a /N)\g) = j a \g), j a being the corresponding eigenvalue. In other 
words, when obtaining some expectation values of intensive operators such as J a /N or H/N the 
ground state \g) can be pictured as a classical angular momentum with fixed coordinates in the 
three-dimensional space (see Fig. 

This point of view makes it clear why such operators ought to be intensive. Otherwise, such 
a classical limit is not valid and terms of order 1 would be important for the calculations of the 
properties of the LMG model. Obviously, all these concepts can be extended to more complicated 
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Hamiltonians such as the extended LMG model, or even Hamiltonians including interactions of 
higher orders as in 114511 . 
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